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RADIO  FREQUENCY  LINAC  DRIVEN  FREE-ELECTRON 
LASER  CONFIGURATIONS 


Introduction 

Significant  advances  in  radio  frequency  (RF)  linac  technology  concerning  the  reduction 
of  emittance  and  increase  in  beam  current  have  been  taking  place.  Lower  emittance  allows 
for  a  smaller  electron  beam  radius  and  thus  higher  current  densities  and  laser  gain.  Under 
appropriate  conditions,  the  radiation  generated  in  a  free  electron  laser  (FEL)  can  be  focused 
and  overcome  the  natural  tendency  to  diffract  (spread  transversely).  This  optical  guiding 
phenomenon  was  first  analyized  in  the  low  gain  regime1,  and  later  in  the  exponential  gain 
regime2-7.  Optical  guiding  in  the  FEL  can  enhance  the  gain  and  efficiency  and  is  an 
integral  part  of  the  configurations  illustrated  in  this  paper. 

We  have  presented  a  formalism  for  describing  the  three-dimensional  radiation  field 
in  FELs  called  the  source  dependent  expansion  (SDE)  method6-7.  In  this  approach,  the 
radiation  field  is  decomposed  into  a  sum  of  complete  normal  modes.  Instead  of  using 
the  usual  modal  expansion  consisting  of  vacuum  Gaussiar  T_  ag  lerre  functions,  the  SDE 
method  incorporates  the  source  function  (driving  current)  s  r  insistently  into  the  func¬ 
tional  dependence  of,  i)  the  radiation  waist,  ii)  the  radiation  wavefront  curvature  and  iii) 
the  complex  radiation  amplitude.  The  SDE  method  significantly  reduces  the  number  of 
modes  that  are  required  to  model  the  radiation  field  and  reduces  the  computational  time 
and  cost  of  the  simulations. 

In  our  numerical  examples  the  electron  dynamics  include:  i)  betatron  oscillations,  ii) 
finite  emittance  and  iii)  energy  spread.  The  FEL  equations  together  with  the  SDE  method 
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are  implemented  in  the  computer  code  SHERA.  The  code  is  applied  to  the  analysis  of  high 
extraction  FELs  driven  by  RF  linacs.  Two  FEL  configurations,  which  we  considered, 
are  the  i)  master  oscillator  power  amplifier  (MOPA)  and  ii)  single  stage  power  amplifier 
(SSPA).  The  schematic  configurations  are  shown  in  Fig.  1.  In  both  of  those  configurations, 
optical  guiding  plays  a  central  role. 

The  Formulation  of  3-D  FEL  Equations 


We  will  consider  a  realistic,  linearly  polarized  wiggler  expressed  in  terms  of  the  vector 
potential, 

Aw(x,y,z)  =  Ax(x,y,z)cos(  /  kw(z')dz')ex,  (1) 

Jo 

where  Ax(x,  y,  z)  and  fcu>(z)  are  slowly  varying  functions  of  z.  The  expressions  for  Az(x,y,  r) 
takes  on  different  forms  depending  on  the  shape  of  the  magnetic  pole  face. 

The  linearly  polarized  radiation  field  is 


A R(r,0,z,t)  =  —  A(r,0,  z)exp[iu>(z/c  -  t)\ex/ 2  +  c.c., 


(2) 


where  A(r,0,z)  —  \A(r,9,  z)\e"i>l'r'0,z^  is  the  complex  amplitude  of  the  radiation  field  ex¬ 
pressed  in  polar  variables,  u>  is  the  frequency,  and  c.c.  denotes  the  complex  conjugate. 
The  wave  equation  governing  A.r  is 

Add  Id2  d2  id2 


A  _  T  A 

r  dr'  dr  '  r2  d92  '  dz 2  c2  dz2  )  R  c  ***  ’ 


r  d-  #.  ft/n 


(3) 


where  Jx  is  the  driving  current  density.  The  detailed  representaion  for  the  current  is  given 
in  Ref.  8.  Using  the  SDE  method6-7,  the  radiation  field  is  represented  in  the  form, 


A(r,^,z)-^Cm>p(^,z)£»P(r,z). 


m,p 


where 


\/2m0c2 


C£(0,z)  =  am>r(z)  cos(p0)  +  bmiP(z)sin{p9), 


(4a) 


m 


Wf)=  (*<•> 


2 


and 


*  = 


2r2 


i; 


r»(z) ' 

=5m,P- 


(4rf) 


(4e) 


The  functional  dependence  of  the  spot  size  r3  and  curvature  a  are  governed  by  the  driving 
current. 

Multiplying  both  sides  of  Eq.  (3)  by  exp[t(u>(z/c  —  f)]  and  integrating  over  t,  from  0 
to  27t/o>,  then  multiplying  both  sides  by  D £,  cos(p0)  or  D ^  sin(p0)  and  integrating  over  r, 
from  0  to  oo,  yields  the  following  SDE  equations  governing  the  radiation  field, 


dz 


4-  Qm,p(Z) 


am,p{Z ) 
bm,p(Z ) 


-  ty  m(m  +  p)B(z) 


—i  A  (m  +  l)(m  +  p  +  1  )B*{z) 


Om+l,p(2) 

&m+l,p(2) 


flm-l,p(2) 

bm-l,p{Z) 

Fm,p(Z ) 

Gm,p(z) 


where 


Qm,P(z )  =  ~  +  i(2m  +  p  +  1) 
r. 


B(z)  =  - 


r's  ,  c  (1  —  a2) 

a - 1- - z - 
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cos  pO 
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\  iflz 
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where  k  =  (1/2)(u>2/c2)(c/ci;)Fb<7(,,  on  top  of  a  quantity  denotes  the  instantaneous 
values  of  the  variable  of  a  given  particle  at  the  axial  position  z,  *  denotes  the  complex 
conjugate  and  the  prime  denotes  a  derivative  with  respect  to  z,  i.e.,  '  =  d/dz.  In  addition, 
aw  =  ( |e  |  /  \Z2m0c2)Ax  is  the  normalized  wiggler  amplitude,  Fb  =  Jo{b)  —  J\{b)  for  a  linearly 
polarized  wiggler,  b  —  au,(0,0,  z)2/2(l4-au,(0,0,  z)2  ),  cr;,  is  the  effective  area  associated  with 
the  transverse  electron  beam  profile,  u >«,  =  (47rje|2no/moc2)1^2  is  the  plasma  frequency  of 
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the  election  beam  on  axis,  6Pt o  =  1  for  p  =  0,  and  6Pto  =  0,  otherwise.  The  symbol  () 
denotes  the  average 


((•••)) 


J  / dx°  J dyo  j dpx0  J dpy0  J dp*°f°Wo'Xo'y°'P*o'Pv°'Pto)(—)'  (9) 


and  is  the  integral  over  all  the  various  initial  conditions  at  z  =  0,  where  f0  is  the  initial 
electron  beam  distribution  function  and  f0  is  normalized  such  that  <  (1)  >  —  1 .  For 
<i  parabolic  transfers'*  electron  beam  profile,  <T(,  =  7r(rj;/2),  and  k  =  u(c /uj)Fb,  where 
v  —  u/l(Tf,/4nc2  ~  I(A)/17  x  103  is  Budker’s  constant  and  rj,  is  the  edge  radius  of  the 
electron  beam. 

The  radiation  amplitude  evolves  as  electrons  move  in  the  longitudinal  ponderomotive 
potential  field  formed  by  the  beating  of  the  wiggler  field  and  the  radiation  field.  The 
electron’s  phase  in  the  ponderomotive  potential  field  is  defined  as  ip,  where 


il'{z-,ip0,xo,y0,px0,p  0,pz0)=  /  (  -  +  kw(z')  -  )dz'  +  4>  +  rp0 , 

Jo  \C  Vz{=')J 


ip o  is  the  initial  phase  in  the  ponderomotive  potential  well,  (xo>2/o)  are  the  initial  trans¬ 
verse  positions,  (pxOiPyo)  are  the  initial  transverse  momentum  and  pz0  is  the  initial  axial 
momentum. 

At  this  point  in  the  SDE  analysis,  the  function  H(z)  is  arbitrary.  If  H(z )  is  not 
specified,  the  equations  for  amiP  and  6miP  in  (5)  are  underdetermined,  i.e.,  there  are  more 
equations  than  unknowns.  The  SDE  method  provides  a  representation  for  r ,  and  a  that 
yields  analytical  results  and  minimizes  the  numerical  computation  in  a  fully  self-consistent 
simulation.  It  is  assumed  that  the  radiation  beam  profile  remains  approximately  Gaussian 
with  a  nearly  circular  cross  section.  In  this  case  we  expect  the  magnitude  of  the  coefficients, 
am  p(r)  and  6m  p(z)  to  become  progressively  smaller  as  m  and  p  take  on  larger  values. 
Hence,  if  the  lowest  order  mode  a o,o  gives  a  rough  approximation  to  the  radiation  field,  we 
may  solve  for  a0io(:),  and  a(z).  We  find 


( djdz  +  C?o,o)ao,o  —  ~iFo,o, 


(ID 


and 


B  —  F],o/a0,0. 


(12) 
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Equation  (12)  above  yields  the  following  first  order  coupled  differential  equation  for  ra{z) 
and  a(z), 

r,.“2“  =  “7f_(Fl’o)/’  (13a) 

lo  r,  a.0,0 


_2J^72  ~  ~  ~  -a(-fl.o)/]  , 

u/  rg  1*0,0 


(136) 


where  ()h,/  denotes  the  real  and  imaginary  parts  of  the  enclosed  function.  In  the  SDE 
method,  the  function  H(z)  is  specified  by  (13). 

The  particle  motion  in  the  z-direction  is  best  written  in  terms  of  the  equation  for  the 
phase  rp  and  the  total  relativistic  gamma, 


and 


where  f3z  —  vz/c. 


d*y  u ;  /  c 

_ - —FBdw(x,y,z)\a(x,y,z)\sin  ip, 

az  7 


^  +  Mt)> 


dz 


/M1  +  /W 


(14) 


(15) 


Description  of  the  3-D  FEL  Code  SHERA 

SHERA  is  a  steady-state,  single  frequency  code,  which  uses  the  SDE  method  to  eval¬ 
uate  the  radiation  beam  in  the  FEL.  It  incorporates  realistic  wiggler  effects,  such  as  trans¬ 
verse  gradient  and  wiggler  noise.  In  the  results  presented  here,  the  wiggler  is  assumed  to 
be  free  of  field  errors  and  it  has  weak  focusing  in  both  the  x  and  y  planes  via  parabolic 
pole  peices9,  i.e., 

Ax(x,y,z )  =  Aw(z)cosh(kw<xx)cosh(kw<yy),  (16) 

where  k^ tX  -f  y  =  k^.  For  equal  focusing  in  the  x  and  y  plane  implies,  kWtX  —  kWtV. 

The  wiggler  field  can  be  tapered  by  assigning  a  functional  form  for  fcu,(z),  aw(z)  or 
resonant  phase  associated  with  a  test  particle.  The  resonant  phase  is  defined  as  the  phase 
associated  with  the  bottom  of  the  ponderomotive  potential  well,  V’fl  =  constant.  A  test 
particle  trapped  at  the  resonant  phase  satisfies  the  following  equation 

~TT'  —  0  =  T(z)  -  R]  sinV’K*  (17) 
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where 

T(s  _  rffcto  _  u/cda2w 
{Z>  dz  2y2  dz  ’ 

and  A';  =  ((1  +  a^)/74)(u>/c)  Fs  Q,w  ja|  is  the  synchrotron  wave  number,  where  de¬ 
notes  the  value  associated  with  the  test  particle. 

For  simplicity,  the  electron  distribution  function  is  assumed  to  be  separable  in  this 
paper,  i.e.,  f0  =  fp{‘<l>o)fe(xo,yo,Pzo,Pyo)fz(Pzo)-  The  initial  transverse  electron  beam  dis¬ 
tribution  is  assumed  to  have  a  water  bag  distribution,  i.e.,  the  electron  density  is 


/<  = 


Rb  <  rb, 
Rb  >  rbi 


(18) 


where  Rl  =  x%  +  yl  +  (v20  +  v20)/(kp(0)c)2  is  the  radius  in  the  emittance  space,  V  is  the 
4-dimensional  phase  space  volume  and  kp  —  aw(x,y,  z)kw /-y  is  the  betatron  wave  number 
for  a  parabolic  pole  face  with  equal  focusing  in  both  the  x  and  y  planes.  Integrated  over 
the  transverse  velocity,  the  electron  beam  has  a  parabolic  density  profile  with  the  beam 
density  n0  on  axis  where  ^  is  the  beam  edge  radius. 

The  transverse  motion  of  the  electrons  due  to  finite  emittance  in  the  realistic  wiggler 
can  be  expressed  accurately  by  an  analytical  expression  for  ykw  <<  1.  The  particle  motion 
is  approximately, 

kpjO) 

b(z) 

where  4*  =  fQz  kp(z')dz'. 

In  the  following  simulations,  the  electron  beam  is  matched  at  the  entrance  of  the 
wiggler,  i.e.,  the  electron  beam  radius  inside  the  wiggler  is  uniform  if  the  wiggler  parameters 
are  constant.  The  matching  condition  becomes 


1/2 


x0  cos  $  +  v*o  sin  $ 
y0  cos  4*  4-  vvo  sin  4> 


(19) 


rb  = 


(20) 


where  Au,  is  the  wavelength  of  the  wiggler,  e„  is  the  normalized  beam  edge  emittance, 
7ren//?7  is  the  electron  beam  area  in  transverse  phase  space. 

The  beam  current  is  modeled  by  discrete  macro  particles.  To  minimize  the  num¬ 
ber  of  electrons  in  the  simulation,  the  integrals  over  the  initial  conditions  are  performed 
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by  quadratures.  Gaussian  quadrature  is  used  for  the  initial  phase  V’o  and  a  72  point 
quadrature10  is  used  for  the  4-dimensional  transverse  emittance.  The  energy  spread  is 
modeled  by  Gaussian,  top-hat  or  any  prescribed  distribution  functions,  and  implemented 
by  either  evenly  spaced  grids,  or  Gaussian  quadrature.  For  simulations  in  this  paper,  the 
full  interval  of  the  Gaussian  quadrature  for  the  energy  spread  is  1.25  times  the  Gaussian 
width. 


RF  Linac  Driven  Power  Amplifiers 

RF  linacs  have  only  been  used  as  drivers  for  FEL  oscillators  due  to  low  gain.  Because  of 
significant  continuous  advances  in  reducing  beam  emittance  and  increasing  beam  current, 
i.e.,  high  beam  brightness,  RF  linacs  appear  to  be  suitable  to  drive  FEL  amplifiers.  With 
the  improvement  in  electron  beam  brightness,  the  high  power  FEL  oscillator  configuration 
becomes  less  attractive,  mainly  because  of  the  problem  of  mirror  damage.  In  this  paper 
we  present  two  configurations  which  can  achieve  high  gain  and  high  efficiency.  These  are 
the  master  oscillator  power  amplifier  (MOPA)  and  the  single  stage  power  amplifier  (SSPA) 
configurations. 


I.  Master  Oscillator  Power  Amplifier  (MOPA) 

The  MOPA  concept,  Fig.  la,  consists  of  trapping  the  electrons  in  the  ponderomotive 
potential  well  and  tapering  the  wiggler  from  the  entrance  to  the  exit  of  the  wiggler  for 
efficiency  enhancement.  To  create  a  large  ponderomotive  potential  well  requires  high  power 
input  radiation  from  a  master  oscillator.  The  large  input  laser  can  be  an  FEL  oscillator. 
We  will  give  an  example  of  an  RF  linac  that  will  produce  1  pm  radiation  in  the  MOPA 
configuration  with  high  efficiency  and  high  power.  The  parameters  of  the  wiggler,  the 
input  radiation  and  the  electron  beam  are  given  in  Table  I. 

The  electron  beam  has  an  energy  of  140  MeV ,  a  peak  current  of  500  A,  an  intrinsic 
energy  spread  of  1.0%  FWHM,  a  normalized  edge  emittance  of  e„  =  0.008  cm  —  rad  and  a 
beam  brightness  of  B  =  1.3  x  1010  A/(m  —  rad)2 .  These  parameters  are  not  too  different 
from  the  existing  RF  linac  beam  at  Boeing  Aerospace  Company. 
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The  amplitude  of  the  vector  potential  of  the  wiggler  is  tapered  for  efficiency  enhance¬ 
ment,  see  Fig.  2.  The  taper  is  slow  initially  and  becomes  faster  at  the  end  of  the  wiggler. 
If  all  the  electrons  are  trapped,  the  maximum  efficiency  would  be  15.8%. 

The  following  simulation  results  are  obtained  with  an  injected  radiation  that  is  focused 
at  the  entrance  of  the  wiggler  with  a  flat  wavefront.  The  efficiency  and  radiation  guiding 
can  be  further  improved  by  moving  the  focal  point,  but  will  not  be  presented  here. 

For  an  input  radiation  power  of  50  MW,  focused  to  a  spot  size  of  0.1  cm,  the  fractional 
trapping  potential  is 


The  estimate  of  the  trapping  protential  has  to  be  adjusted  by  the  initial  resonant  phase, 
which  is  sinV'fi  =  0.25.  The  actual  trapping  potential  is  sufficient  for  trapping  an  electron 
beam  with  a  1.0%  energy  spread.  Since  the  radiation  grows  rapidly,  some  of  the  electrons 
not  originally  inside  the  separatrix  can  become  recaptured. 

In  this  illustration,  an  axially  symmetric  electron  beam  and  radiation  beam  are  as¬ 
sumed.  The  wiggler  has  parabolic  pole  faces  so  that  weak  transverse  focusing  is  the  same 
in  both  the  x  and  y  planes.  A  total  of  10  Gaussian-Laguerre  modes  were  used  in  the  sim¬ 
ulation.  The  differences  in  results  using  10  Gaussian-Laguerre  modes  versus  6  Gaussian- 
Laguerre  modes  is  typically  much  less  than  5%. 

Figure  3  is  a  plot  of  the  radiation  spot  size  as  a  function  of  distance.  In  the  absence 
of  optical  guiding,  the  Rayleigh  length  would  be  3  in,  and  at  the  end  of  the  interaction 
the  radiation  spot  size  would  be  13  times  larger  than  the  initial  spot  size  of  0.1  cm.  As  a 
result  of  optical  guiding,  the  final  spot  size  is  only  0.32  cm.  The  radiation  is  therefore  able 
to  propagate  through  the  vacuum  chamber  without  significant  loss,  because  the  condition 
g  —  2.5  cm  >  5 r,  +  0.3  cm  =  1.9  cm  is  satisfied,  where  g  is  minimum  wiggler  gap11. 

Figure  4  is  a  plot  of  the  radiation  power  in  each  mode  as  a  function  of  axial  distance. 
The  fundamental  mode  clearly  dominates  and  the  SDE  method  is  extremely  accurate  for 
evaluating  the  spot  size  and  wave  front  curvature. 

Figure  5  is  a  plot  of  the  efficiency  as  a  function  of  beam  energy  spread.  For  electron 
beams  with  energy  spread  less  than  1.0%  Gaussian  full  width,  the  smaller  input  radiation 
power  gives  a  higher  efficiency.  This  is  a  consequence  of  optical  guiding.  We  have  shown 
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in  Ref.  6  that  perfectly  guided  radiation  beams  cannot  be  maintained  in  the  high-gain 
trapped-particle  regime.  However,  the  radiation  is  well-focused  compared  to  free  space 
diffraction.  The  envelope  equation  for  the  radiation  beam  has  the  form6-7 

r'J  +  K2ra  =  0.  (22) 

The  parameter  governing  radiation  focusing  is 

2 

A'2  =  ^ ^  ^  —  1  4-  C2  <  sint/>  >2  4-2 C  <  cosip  >  -f—  r2C'  <  simp  ,  (23) 

where  C(z)  =  —  /)(1  4-  /)_7(al„/|ao,o(2:)|)  is  the  coupling  coefficient,  f(z)  = 

cr\y/<Tr  is  the  filling  factor,  and  <rr  =  irr2  is  the  area  of  the  radiation  beam.  The  first 
term  on  the  right-hand  side  of  (23)  is  defocusing  and  corresponds  to  the  usual  diffraction 
expansion,  the  second  and  third  terms  are  always  focusing  while  the  last  term  is  usually  a 
defocusing  contribution.  For  perfectly  guided  beams,  K 2  =  0.  For  weakly  guided  beams, 
the  focusing  terms  are  reduced.  The  coupling  coefficient  C  is  larger  when  the  initial  input 
is  smaller.  Thus,  tue  radiation  focuses  better  when  the  input  power  is  smaller.  Figure  6 
is  a  comparison  of  radiation  spot  sizes  as  a  function  of  z  for  the  input  radiations  of  ( — ) 

50  MW  and  ( - )  bOOMW  obtained  for  simulations  with  no  energy  spread.  The  spot  size 

associated  with  50  MW  input  is  smaller  and  the  final  efficiency  is  larger. 

II.  Single  Stage  Power  Amplifier  (SSPA) 

A  somewhat  simpler  configuration,  that  can  achieve  high  power  and  high  efficiency, 
is  the  single  stage  power  amplifier  (SSPA)12.  In  the  SSPA  configuration  (see  Fig.  lb),  a 
low  power  laser,  such  as  a  dye  laser,  is  directly  amplified  and  optically  guided  in  a  uniform 
wiggler  region  until  saturation  occurs.  The  intrinsic  efficiency  in  the  uniform  wiggler 
region  is  small,  typically  less  than  2%.  For  efficiency  enhancement,  the  wiggler  is  spatially 
tapered  just  prior  to  saturation.  At  this  point  the  electron  beam  is  tightly  bunched  in 
the  ponderomotive  wave,  and  the  trapping  potential  can  be  sufficiently  large.  The  SSPA 
configuration  requires  a  slightly  longer  overall  wiggler  length.  The  main  advantage  of  the 
SSPA  configuration  is  the  elimination  of  the  high  power  master  oscillator. 
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The  parameters  of  the  wiggler,  the  input  radiation  and  the  electron  beam  for  the  1  pm 
SSPA  configuration  are  given  in  Table  II.  The  electron  beam  has  an  energy  of  250  MeV , 
a  peak  current  of  1  kA,  an  intrinsic  energy  spread  of  1.0%  FWHM,  a  normalized  edge 
emittance  of  e„  =  0.007  cm  rad  and  thus,  a  beam  brightness  B  =  5  x  1011  A/(m  -  rad)2 . 

The  following  simulation  results  are  obtained  with  an  injected  radiation  beam  power 
of  15  kW  peak,  focused  at  the  entrance  of  the  wiggler  to  a  spot  size  of  0.11  cm.  The  free 
space  Rayleigh  length  is  4  m.  At  exact  resonance,  the  radiation  saturates  at  a  distance  of 
15  m  with  an  intrinsic  efficiency  of  0.39%.  At  the  end  of  the  uniform  wiggler  region,  the 
trapping  potential  is 

\e\ttrap  =  4  /  W1/2  =  3  2% 

7 m0c2  V  1  +  a2,  ) 

In  the  uniform  wiggler  section,  the  radiation  beam  is  optically  guided  and  the  magnetic 
field  amplitude  and  period  are  3.155  kG  and  9.6  cm,  respectively.  The  wiggler  taper  begins 
just  prior  to  the  saturation  point  (z  =  15  m).  Figure  7  shows  the  taper  of  the  amplitude  of 
the  vector  potential  of  the  wiggler.  If  all  the  electrons  are  trapped,  the  maximum  efficiency 
would  be  27%.  A  plot  of  the  efficiency  as  a  function  of  the  axial  distance  is  shown  in  Fig. 
8.  The  maximum  efficiency  is  22%  at  the  end  of  the  wiggler  and  a  substantial  fraction 
of  the  beam  electrons  are  trapped  in  the  ponderomotive  buckets.  Figure  9  indicates  that 
the  wiggler  can  be  tapered  further  to  obtain  even  higher  efficiencies.  The  spot  size  of  the 
radiation  field  is  plotted  in  Fig.  9.  The  spot  size  at  the  end  of  the  wiggler  is  0.24  cm. 

Conclusions 

With  the  recent  advances  in  RF  linac  technology  towards  higher  power  and  higher 
brightness,  FEL  configurations  driven  by  RF  linacs  are  no  longer  limited  to  oscillators. 
We  have  presented  exan.r  -s  of  the  MOPA  and  SSPA  configurations  for  achieving  high 
FEL  power  and  efficient .  xhe  SSPA  is  the  simpler  of  the  two  configurations  because  the 
high  power  master  oscillate-  ieplaced  with  a  low  power  conventional  laser  source. 
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Table  I 


Parameters  of  Numerical  Simulation  in  the  MOPA  Configuration 


Electron  Beam  (parabolic  density  profile) 
Energy,  E 
Current ,  I 

Intrinsic  E~spread,  AE/E 
Slippage  distance 
Norm,  edge  emittance,  en 
Radius ,  r 

Energy  spread  due  to  en, 

A-r,  _  1 


=  r* 


2  1 +al 


Wiggler  (linearly  polarized) 

Period,  \w 

Init.  magnetic  field,  Bw 
Final  magnetic  field,  Bw 
Wiggler  length,  Lw 

Initial  normalized  vector  potential,  aw 
Taper  (amplitude) 


140  MeV 
500  A 

1.0%  (Gaussian  full  width) 
0.08  cm 
0.008  cm-rad 
0.08  cm 

0.17% 


5.0  cm 
4.28  kG 
3.21  kG 
40  m 

V2 

non-uniform  (max  15.8%) 


Radiation 

Wavelength,  A 

Input  power ,  P\n 

Min.  spot  size  of  input  rad. 

Initial  Rayleigh  length 

Final  spot  size,  r,  (z=40  m) 

Output  power,  Pout  (peak) 

Efficiency  at  end  of  wiggler 


1  fim 
50  MW 

0.1  cm  (at  z  =  0  m) 

3.0  m 

0.32  cm 

5.2  GW 

7.5% 
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Table  II 


Parameters  of  Numerical  Simulation  in  the  Power  Amplifier  Configuration 


Electron  Beam  (parabolic  density  profile) 
Energy,  E 
Current ,  / 

Intrinsic  E-spread,  A  E/E 
Slippage  distance 
Norm,  edge  emittance,  tn 
Radius ,  r<, 

Energy  spread  due  to  en, 

Alt  -  Ill  ’ 

7,  •> 


2  r l  1+aJ, 


Wiggler  (linearly  polarized) 

Period,  A^ 

Init.  magnetic  field,  Bw 
Final  magnetic  field,  Bw 
Wiggler  length,  Lw 

Initial  normalized  vector  potential,  aw 
Taper  (amplitude) 


250  MeV 
1000  A 

1.0%  (Gaussian  full  width) 
0.08  cm 
0.007  cm-rad 
0.087  cm 

0.06% 


9.6  cm 
3.155  kG 
2.05  kG 
80  m 
2.0 

non-uniform  (max  27%) 


Radiation 

Wavelength,  A 

Input  power,  Pln 

Min.  spot  size  of  input  rad. 

Initial  Rayleigh  length 

Final  spot  size,  r,  (z=80  m) 

Output  power,  Pout  (peak) 

Efficiency  at  end  of  wiggler 


1  fxm 

15  kW 

0.11  cm  (at  z  =  0  m) 
4.0  m 
0.23  cm 
53  GW 
22% 
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Fig.  3.  Plot  of  the  radiation  spot  size  as  a  function  of  distance  z  for  1.0%  Gaussian 
full  width  energy  spread. 
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Fig.  4.  Plot  of  radiation  power  for  each  mode  as  a  function  of  distance  with  initial 
radiation  input  power  of  50  MW. 
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Fig.  5.  Plot  of  the  Efficiency  as  a  function  of  energy  spread  for  input  radiation  power 
( — )  50  MW  and  ( - )  500  MW. 
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Fig.  6.  Plots  of  spot  sizes  as  a  function  of  axial  distance  for  an  electron  beam  with 
no  initial  energy  spread  for  input  radiation  power  of  ( — )  50  MW  and  ( - )  500  MW. 
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NORMALIZED  VECTOR  POTENTIAL 
OF  THE  WIGGLER,  ^ 


Fig.  7.  Plot  of  the  amplitude  of  the  normalized  vector  potential  of  the  wiggler  as  a 
function  of  interaction  distance  z. 
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Fig.  8.  Plot  of  the  radiation  spot  size  as  a  function  of  distance  z  for  1.0%  Gaussian 
full  width  energy  spread. 
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Fig.  9.  Plot  of  efficiency  as  a  function  of  the  axial  distance  z. 
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